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Ice- Affected  Flood  Frequency 


Virtually  every  hydraulic  structure 
design  or  other  type  of  study  involv¬ 
ing  rivers  requires  the  development 
of  a  frequency  analysis,  in  which 
river  discharge  is  plotted  against 
exceedance  probability.  These  curves 
allow  the  quantitative  selection  of  de¬ 
sign  life  for  structures  such  as  bridges 
or  levees,  and  are  required  for  risk- 
based  analyses.  When  combined  with 
stage-discharge  and  stage-damage 
curves,  damage  estimates  and  other 
economic  analyses  can  be  made. 
Stage-frequency  curves  also  are  often 
the  basis  for  regulatory  floodplain 
delineations  such  as  those  presented 
in  flood  insurance  studies. 

Guidance  is  available  for  the  prep¬ 
aration  of  discharge-frequency  curves 
for  open-water  conditions  at  gaged 
sites  (e.g..  Water  Resources  Council 
1982,  Hydrologic  Engineering  Center 
1992). 

Because  ice  jams  often  cause  high 
stages  at  relatively  low  discharges, 
the  traditional  discharge-frequency 
curve  is  not  appropriate  and  a  stage- 
frequency  curve  must  be  constructed. 
Developing  an  ice-affected  stage-fre¬ 
quency  curve  can  be  difficult.  Gages 
subject  to  ice  jams  may  experience 
damage  that  renders  them  useless 

(Fig- 1). 

Ice  jams  are  unstable  at  high  dis¬ 
charges,  and  therefore  may  affect  only 
a  portion  of  the  stage-frequency 
curve.  Some  stages  reported  to  be 
ice-affected  may  actually  be  open-wa¬ 
ter  stages.  Length  of  record,  the  ma¬ 
jor  obstacle  encountered  in  open-wa¬ 
ter  flood  frequency  analyses  at  gaged 
sites  (Greis  1983),  is  exacerbated  for 
ice-related  events,  which  generally 
occur  much  less  frequently  than  open- 


water  events  and  consequently  have 
even  smaller  sample  sizes. 

A  recent  study  by  the  U.S.  Army 
Corps  of  Engineers  St.  Paul  District 
(USAGE  1999a,b)  provides  a  good 
example  of  developing  ice-affected 
stage  frequency  along  a  river  reach, 
hollowing  the  large  flood  of  1997,  it 
was  necessary  to  update  the  flood 
insurance  studies  for  the  communi¬ 
ties  of  Wahpeton,  North  Dakota,  and 
Breckenridge,  Minnesota,  to  include 
the  effects  of  ice,  upstream  flood 
control  reservoirs,  and  an  upstream 
breakout-flow  area  (i.e.,  above  cer¬ 
tain  elevations,  overland  flow  will 
b5q)ass  the  gage). 

The  analysis  presented  here 
addresses  the  development  of  ice- 
affected  stage  frequency.  The  method¬ 
ology  employed  is  within  the  general 


guidance  for  ice-affected  flooding  as 
provided  by  the  Eederal  Emer¬ 
gency  Management  Agency  (EEMA 
1995).  Specific  guidance  was  pro¬ 
vided  by  the  Corps  of  Engineers 
Cold  Regions  Research  and  Engi¬ 
neering  Laboratory  (CRREL),  Han¬ 
over,  New  Hampshire,  and  the  Corps 
of  Engineers  Hydrologic  Engineering 
Center  (HEC),  Davis,  California. 

Site  description 

The  Red  River  of  the  North  is 
formed  by  the  confluence  of  the  Bois 
de  Sioux  River  and  the  Otter  Tail 
River  in  the  communities  of  Wahpe¬ 
ton,  North  Dakota,  and  Brecken¬ 
ridge,  Mirmesota.  The  local  effective 
drainage  area  downstream  of  two 
flood  control  reservoirs  (White  Rock 
Dam  and  Orwell  Dam)  is  1020 


Figure  1.  Gage  house  for  USGS  Gage  No.  01010000,  St.  John  River  at  Ninemile  Bridge, 
Maine,  following  April  1991  ice  jam.  Note  remaining  ice  debris,  right. 
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square  miles.  A  general  location 
map  is  shown  in  Figure  2. 

Gaged  streamflow  data  is  pro¬ 
vided  by  uses  Gage  No.  05051500, 
Red  River  of  the  North  at  Wahpeton. 
The  gage  is  located  800  ft  down¬ 
stream  from  the  confluence  of  the 
Bois  de  Sioux  River  and  the  Otter  Tail 
River  at  River  Mile  548.6.  The  period 


of  record  dates  from  April  1942  with 
observed  flow  measurements  that 
are  considered  fair.  The  maximum 
instantaneous  peak  stage  of  19.42  ft 
was  influenced  by  backwater  from 
ice  and  occurred  on  6  April  1997. 
The  hydrograph  of  fhe  1997  even!  is 
shown  on  the  stage-discharge  curve 
at  the  gage  (Fig.  3).  The  rising  limb  of 


the  hydrograph  generally  follows 
the  ice-affected  armual  peaks  until 
some  point  above  about  8000  cfs, 
after  which  it  rises  to  the  peak  stage 
of  962.5  ft  along  an  apparent  open- 
water  curve.  The  falling  limb  fol¬ 
lows  the  same  open-water  curve. 

Stage  frequency  at  gage 

Because  of  fhe  significanf  number 
of  ice-affeefed  floods  in  fhe  study 
area,  flood  frequency  analysis  based 
on  peak  annual  instantaneous  dis¬ 
charges  alone  was  not  considered 
appropriate.  Of  the  56  years  of  avail¬ 
able  sfreamflow  records  (1942-1997) 
for  the  Red  River  of  the  North  at 
Wahpeton,  half  of  fhe  annual  peak 
stages  were  ice-affected  and  half 
were  associafed  wifh  open-water  dis¬ 
charges.  Using  the  mixed-population 
approach,  an  accepted  method  of 
analyzing  sfages  arising  from  differ¬ 
ent  causes,  the  two  populations  were 
separated  according  to  the  season 
(ice-affected  vs.  open-water)  and 
not  arbitrary  calendar  months,  as 
is  sometimes  done.  In  the  mixed- 
population  method,  the  annual  peak 
ice-affected  stage  frequency  and 
annual  peak  open-water  stage 
frequency  are  combined  using 

Pc  =  Pl+  Pl-iPdiPl)  (1) 
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Figure  3.  Elevation-discharge  data  for  USGS  Gage  05051500,  Red  River  of  the  North  at 
Wahpeton,  North  Dakota,  showing  ice-affected  and  annual  peaks  along  with  data  for  the 
1997  event  and  the  current  rating  curve  (#17). 


where  P,.  is  fhe  probability  of  a 
selecfed  elevation  being  equaled  or 
exceeded  from  either  an  ice-affected 
flood  even!  or  an  open-water  flood 
evenf,  is  the  probability  of  fhaf 
elevafion  being  equaled  or  exceeded 
from  an  ice-affected  flood  event,  and 
P2  is  the  probability  of  fhaf  elevafion 
being  equaled  or  exceeded  from  an 
open-water  flood  event. 

The  mixed  population  method 
requires  that  peak  stages  be  deter¬ 
mined  for  each  season  for  every  year 
of  record.  It  can  be  difficult  to  locate 
such  information  for  ice-affeefed 
gages,  and  ice  peaks  for  years  in 
which  fhe  annual  peak  was  not  ice- 
affected  must  often  be  synthesized. 
Fortunately,  in  this  case  a  detailed 
search  of  archived  USGS  records 
revealed  ice  peaks  for  fhe  years  in 
which  fhe  annual  peak  occurred  dur¬ 
ing  the  open  water  season.  The  ice- 
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Table  1.  Stage  frequency  for  the  Red  River  of 
the  North  at  Wahpeton. 

Elevation 

Exceedence  probability 

(ft) 

Ice-affected 

Open-water 

Combined 

964.8 

0.002‘ 

964.1 

0.003‘ 

963.0 

0.0045 

0.0025 

0.0070 

962.0 

0.0080 

0.0050 

0.0130 

961.0 

0.0150 

0.0100 

0.0250 

960.0 

0.0250 

0.0170 

0.0416 

959.0 

0.0430 

0.0300 

0.0717 

958.0 

0.0700 

0.0550 

0.1212 

‘Elevations  for  the  0.2%  and  0.3%  exceedence  frequency 
events  were  defined  by  the  all-season  peak  discharge- 
frequency  relationship  using  the  open-water  elevation- 
discharge  rating  curve. 


affected  stage-frequency  curve  was 
developed  at  the  USGS  gage  based 
on  the  56  years  of  stage  data  avail¬ 
able.  Annual  peak  ice-affected  eleva¬ 
tions  were  graphically  plotted  using 
median  plotting  positions.  The  1997 
peak  ice-affected  elevation  of  962.4  ft 
was  plotted  as  the  highest  flood 
elevation  since  the  historic  spring 
flood  of  1897  (peak  elevation  960  ft) 
and  has  a  median  plotting  position  of 
0.69%  based  on  a  101-year  historic 
period  from  1897  to  1997. 

The  open-water  stage-frequency 
curve  was  based  on  an  annual  series 
peak  open-water  discharge-fre¬ 
quency  curve  over  a  56-year  period. 
An  open-water  elevation-discharge 
rating  curve  developed  at  the  USGS 
gage  was  used  in  combination  with 
the  adopted  open-water  discharge- 


frequency  curve  to  graphically 
develop  an  annual  peak  open- 
water  stage-frequency  curve. 

The  combined-population 
stage-frequency  curve  was 
developed  at  the  gage  by  com¬ 
bining  the  independent  fre¬ 
quency  curves  for  the  ice- 
affected  elevations  and 
open-water  elevations  using 
eq  1.  Because  of  expected 
ice  instability  at  high  dis¬ 
charges,  the  upper  end  of  the 
frequency  curve  (0.2%  and 
0.3%  exceedence  frequency 
elevations)  was  defined 
from  the  annual  all-season 
peak  discharge-frequency 
curve.  The  adopted  combined-popu¬ 
lation  stage-frequency  curve  for  the 
Red  River  of  the  North  at  the 
USGS  gage  is  shown  in  Figure  4 
along  with  the  frequency  curves  for 
ice-affected  elevations  and  open- 
water  elevations.  Table  1  provides  a 
summary  of  exceedence  probabilities 
for  selected  elevations  for  the  ice- 
affected  population,  open-water 
population,  and  the  combined  popu¬ 
lation  for  the  Red  River  of  the  North 
at  the  USGS  gage. 

Stage  frequency 
along  study  reach 

Stage-frequency  curves  were 
computed  at  cross-section  locations 
for  floods  with  selected  recurrence 
intervals  of  10, 50, 100,  and  500  years. 


Exceedence  Frequency (%) 

Figure  4.  Stage-frequency  curves  for  the  Red  River  of  the  North  at  Wahpeton. 


These  in  turn  were  used  to  establish 
the  100-  and  500-year  floodplain 
boundaries  and  the  100-year  flood¬ 
way  throughout  the  study  area. 

The  HEG-2  step-backwater  com¬ 
puter  program  (USAGE  1990)  was 
used  to  generate  water  surface  eleva¬ 
tions  for  selected  recurrence  intervals 
for  both  open-water  conditions  and 
ice-affected  conditions  at  all  cross 
sections  throughout  the  study  area. 
The  HEG-2  model  was  calibrated 
using  the  historical  flood  profiles  of 
the  Red  River  of  the  North,  Otter 
Tail,  and  Bois  de  Sioux  Rivers.  Water 
surface  profiles  at  selected  recur¬ 
rence  intervals  for  open-water  condi¬ 
tions  were  computed  using  the  start¬ 
ing  water  surface  elevations  from 
the  stage-discharge  rating  curve. 

Water  surface  profiles  were  com¬ 
puted  for  ice-affected  conditions 
with  the  USAGE  computer  program 
NEWTHK,  a  version  of  IGETHK 
(Tuthill  et  al.  1998)  that  interfaces 
with  HEG-2  in  an  iterative  manner. 
NEWTHK,  which  has  algorithms  for 
ice  jam  stability  analysis  (Tuthill  and 
Mamone  1998),  allows  the  modeler 
to  determine  the  potential  severity 
of  ice  jamming  to  an  area.  Data 
necessary  for  the  ice  portion  of  the 
hydraulic  analysis  were  obtained 
from  field  observations  at  Wahpeton 
and  Breckenridge  beginning  in  1990. 
The  ice  type  varies  somewhat  from 
year  to  year  but  is  basically  strong 
and  competent.  The  typical  ice  thick¬ 
ness  of  18  in.  was  used  in  the  study. 

A  Manning's  n  of  0.026  for  the  un¬ 
derside  of  the  ice  cover  was  as¬ 
sumed.  The  discharges  used  in  de¬ 
veloping  the  ice-affected  water 
surface  profiles  were  determined 
from  the  ice-affected  stage-discharge 
curve  (Eig.  3)  for  the  stages  on  the 
stage  frequency  associated  with  the 
selected  recurrence  intervals.  Using 
these  discharges,  ice-affected  stages 
at  upstream  reaches  were  deter¬ 
mined  using  NEWTHK. 

The  calibrated  HEG-2  was 
checked  against  the  stage-discharge 
hydrograph  recorded  during  the 
spring  flood  of  1997  (Eig.  4)  to  deter¬ 
mine  how  well  the  calculated  water 
surface  elevations  compared  to 
the  observed  values  at  the  USGS 
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gage.  Over  the  entire  range  of  ice- 
cover  measurements  and  up  to  the 
peak  stage,  on  the  rising  (ice-af¬ 
fected)  and  falling  (open-wafer) 
limbs  of  fhe  hydrograph,  fhe  com- 
pufed  and  observed  values  were 
wifhin  0.2  ff. 

For  each  cross-secfion  locafion, 
frequency  relationships  were  devel¬ 
oped  for  bofh  fhe  ice-affecfed  and 
open-wafer  populafions  by  graphi¬ 
cally  tiffing  curves  fo  fhe  elevafions 
generafed  by  fhe  HEC-2  model  for 
floods  of  eighf  differenf  recurrence 
infervals  (5, 10,  20,  50, 100, 125,  200, 
and  500  years).  The  combined-popu¬ 
lation  sfage-frequency  curve  af  each 
of  fhe  selecfed  cross-secfion  locations 
was  fhen  defermined  using  eq  1. 

Summary 

Engineers  accusfomed  fo  develop¬ 
ing  discharge-frequency  curves  along 
rivers  are  offen  confused  when  con- 
fronfed  wifh  fhe  need  fo  perform  an 
ice-affecfed  frequency  analysis, 
because  sfage  frequency  musf  be  con¬ 
sidered.  This  case  sfudy  provides  an 
example  of  developing  such  infor- 
mafion  using  a  mixed-populafion 
analysis  under  ideal  conditions  (i.e., 
a  USGS  gage  wifh  long-ferm  records 
is  locafed  nearby,  peak  sfages  are 
abouf  equally  divided  befween  open- 
wafer  and  ice-affecfed  sfages,  and 
ice-affecfed  peak  sfages  can  be  iden¬ 
tified  for  each  year  of  record).  These 
are  fhe  essenfial  sfeps  in  such  an 
analysis: 

1.  Defermine  open- wafer  sfage 
frequency  for  all  years  of  record. 

2.  Defermine  ice-affecfed  sfage 
frequency  for  all  years  of  record. 

3.  Combine  dafa  from  sfeps  1  and 
2  for  selecfed  sfages  fo  develop 
an  all-season  sfage-frequency 
curve. 

4.  Perform  open-wafer  backwafer 
analysis  using  discharges  defer¬ 
mined  from  fhe  open-wafer 
sfage-discharge  and  sfage- 
frequency  curves  for  selecfed 
recurrence  infervals  (use  cali- 
brafed  HEC-2  or  HEC-RAS). 

5.  Perform  ice-affecfed  backwafer 
analysis  using  discharges  defer¬ 
mined  from  fhe  ice-affecfed 
sfage-discharge  and  sfage- 


frequency  curves  for  selecfed 
recurrence  infervals  (use  cali- 
brafed  NEWTHK,  ICETHK,  or 
HEC-RAS). 

6.  Combine  fhe  resulfs  of  sfeps 
4  and  5  for  selecfed  sfages  af 
fhe  selecfed  cross  secfions  fo 
develop  an  all-season  sfage- 
frequency  curve. 

The  same  general  approach  may 
be  utilized  under  non-ideal  condi¬ 
tions.  However,  addifional  sfeps  such 
as  synfhesis  of  ice-affecfed  peak 
sfages  for  every  year  of  record  could 
be  required  because  armual  ice  peaks 
are  nof  offen  available. 
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Information  Exchange  Bulletin 

The  Ice  Engineering  Information  Exchange 
Bulletin  is  published  in  accordance  with 
AR  25-30  as  one  of  the  information  ex¬ 
change  functions  of  the  Corps  of  Engi¬ 
neers.  It  is  primarily  intended  to  be  a 
forum  whereby  information  on  ice  engi¬ 
neering  work  done  or  managed  by  Corps 
field  offices  can  be  disseminated  to  other 
Corps  offices,  other  U.S.  Government 
agencies,  and  the  engineering  community 
in  general.  The  purpose  of  the  Ice  Engi¬ 
neering  Information  Exchange  Bulletin  is 
information  exchange  and  not  the  pro¬ 
mulgation  of  Corps  policy;  thus,  guid¬ 
ance  on  recommended  practice  in  any 
given  area  should  be  sought  through 
appropriate  channels  or  in  other  docu¬ 
ments.  This  bulletin's  contents  are  not  to 
be  used  for  advertising,  publication,  or 
promotional  purposes.  Citation  of  trade 
names  does  not  constitute  an  official  en¬ 
dorsement  or  approval  of  the  use  of  such 
commercial  products. 
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Tatinclaux  (CEERD-RI),  72  Lyme  Road, 
Hanover,  NH  03755-1290,  or  call  603- 
i  646-4361.  i 
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